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^ ; Abstract 

Q,^ . A lattice calculation of the ttA^ sigma term is described using dynamical staggered fermions. Preliminary 

results give a sea term comparable in magnitude to the valence term. 

^ \ Introduction and Theoretical Discussion 

^ O ■ Hadrons appear to be far more complicated than the (rather successful) constituent quark model 

\ would suggest. For example a recent result from the EMC experiment, [Q, suggests that constituent 

quarks are responsible for very little if any of the nucleon spin. Another, older result, is from the vrA^ 
sigma term, see e.g. which seems to give a large strange component to the nucleon mass. To explain 
theoretically these results and other experiments involves the computation of certain matrix elements - 
a non-perturbative calculation. This is thus an area where lattice calculations may be of some help. 

The ttA^ sigma term, a^^N , is defined|^ as that part of the mass of the nucleon (for definiteness the 
proton) coming from the vacuum connected expectation value of the up (u) and down (d) quark mass 
terms in the QCD Hamiltonian, 

a^fq ^ 'm{N\uu + dd\N), (1) 

, where we have taken these quarks to have equal current mass {— m). Other contributions to the nucleon 

0^ ' mass come from the chromo-electric and chromo-magnetic gluon pieces and the sea terms due to the 

s quarks. Experimentally this matrix element has been measured from low energy -k-N scattering. A 
delicate extrapolation to the chiral limit || gives a result for the isospin even amplitude of S//^ with 
S = cr-nN, from which the -kN sigma term may be found. The precise value obtained this way has been 
under discussion for many years. However within the limits of our lattice calculation, this will not concern 
, us here and for orientation we shall just quote a range of results from later analyses of a-j^N ~ 56MeV, 

II, down to 45MeV, §. 

To estimate valence and sea contributions to (XTrAf, classical current algebra analyses assume octet 
dominance and make first order perturbation theory about the S'C/i?(3) flavour symmetric Hamiltonian. 
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H ! We have 



iirn = rnluu + dd) + nisSS = rriav^I^ —1= — ipX^-ip- (2) 

v3 

The first term in the second equation is the flavour symmetric Hamiltonian with mass ruav = {'2'ni + 
ms)/V3, while the second term is the flavour breaking piece, {ip is the column vector {u, d, s).) We thus 
have from eq. (|l|) 

CT^JV ~ ■m{N\uU + dd - 2ss\N)syrnm +'2rn{N\ss\N)symm 

= cr^AT + O^ttJV) («>) 

where we have first assumed that the nucleon wavefunction does not change much around the symmetric 
point. We then subtract and add a strange component. At the symmetric point the u and d quarks each 
have equal valence and sea part, while the s quark matrix element only has a sea component. Thus in 
the first term the sea contribution cancels, justifying the definitions given in the last line of eq. (H). Also 
this means that in the following we shall interchangeably talk about 'strange' and 'sea' contributions to 
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^Actually, originally as the matrix element ofjjre double commutator of the Hamiltonian with two axial charges. However 
this is equivalent to the definition given in eq. (hi), see for example ref. 



the nucleon. At first order perturbation theory crJJ'^ easily be calculated* to give cr^"^ « 25MeV and 
so af^ « 31 ~ 20MeV. This in turn means that 

ms{N\ss\N) « 400 - 250MeV, (4) 

which would indicate a sizeable portion of the nucleon mass comes from the strange quark contribution. 
(Remember that the mass of the proton is about 938MeV.) This would not be expected from the 
constituent quark model, where the nucleon is made only of (dressed) u and d quarks. 

Staggered Fermions 

As we see from above, we need to calculate certain matrix elements. One of the most promising 
methods at present to calculate them comes from the lattice, where Euclidean space-time is discretised 
with a lattice spacing a. The problem is then turned into a statistical mechanical one of evaluation of 
correlation functions (as to be described below). This can be attacked using Monte Carlo simulations. 
The continuum limit is obtained when a — > 0, or /3 = 6/.g^(a) — > oo. This procedure is well known, see for 
example Q . We have performed a large scale simulation on a 16^ X 24 lattice of QCD using the Hybrid 
Monte Carlo method. The main emphasis was for (3 = 5.35, m — 0.01, although some other runs were 
also performed. (Our results and further details are given in |^.) 

There are two commonly used formulations of fermions on the lattice - Wilson and staggered. Both 
have their advantages and disadvantages. We have used the staggered approach, which describes in 
the continuum 4 degenerate flavours. Here the fermion flelds in the partition function, x°'{^) ^-^e only a 
function of the lattice point and colour index. Quark spin and flavour combinations are given by summing 
appropriately over x's sitting on the vertices of hypercubes, thus 

qc.f{X)^Y.^tfX{X + A). (5) 

A 

The hypercube is defined hy x ^ X + A, being or 1 and even, where x is the lattice site. 
is given by . . .7^*. Eq. ^ projects out the flavour (/) and spin (a) degrees of freedom. In the 
continuum limit we have for m = a chiral symmetry £'[7^(4) (g) S'?7jj(4). The diagonal symmetry, valid 
for all m, is the flavour group SUf{4:)- The advantage of the staggered lattice fermion formulation is 
that we have, for m = 0, a remnant U (1) (E) U{1) symmetry. This allows a clean study of chiral symmetry 
breaking as m —> 0, and also the emergence of a massless tt, with oc m. (For Wilson fermions, 
one must carefully tunc the hopping parameter to obtain a massless tt, chiral properties always being 
explicitly broken.) As we saw in the last section we have at the symmetric point an SUp{3) symmetry 
for the u, d, s quarks. (Corrections for the strange quarks are treated in first order perturbation theory, 
as part of the sea contribution.) As we shall be here interested in the questions of sigma/nucleon and 
valence/sea ratios, we expect that the different number of quarks present will make little difference. 

Lattice calculations have to be performed when the lattice spacing a is small enough that finite lattice 
effects are negligible. Staggered fermions have a natural indicator to show this by when the fiavour 
symmetry is restored. On the lattice meson {M xx) ^-^d baryon {B XXX) operators which in the 
continuum describe degenerate states lie in several distinct lattice representations, each with a different 
mass. As a — !■ these masses must tend to a common mass. Practically we say that we are close enough 
to the continuum when possible mass differences are negligible. For our configurations we seem to be 
entering such a region, 

Another technical point concerns the existence of a parity partner in the correlation functions. Masses 
{Ma) and amplitudes {Aaa) are measured from 2-point correlation functions. For the baryon we have 

C{t) {B{ t)Bw (0)) 

^Baryon octet mass splittings are given by 

Mb - Msymm = -(rris - m){B\-ip\^ 1p\B) symm / V^, 

where B is the i"*" nucleon multiplct. From the Wigner-Eckart theorem for SU{3) we have 

(B\'^X^-^\B)symm = FtriB^' [X^ , B]) + Dtr(Bt{A8, b}). 

From the known baryon masses _F, D can then be estimated. (We use the numbers given in as {rris — m)F 190, 
{nis — m)D —61.) Furthermore 

ul°^ = m.VZ{N\i,\'^'4)\N)symm = m{3F - D). 
Using ms/m 26 gives the quoted results. 



« {0\BSlBw\0) 

« ANNfiN + Aaa^a iT>i>0. (6) 

In the second line we have gone from the (numerically) calculable 2-point correlation function to the 
equivalent expression in the operator formalism. The transfer matrix is denoted by f = = exp {-2H). 
A shift of one unit still lies in the hypercube, which then involves a flavour transformation. Thus we can 

write iS'4 — S4T2 where S4 is a non-local flavour transformation. (No explicit representation of S4 is at 
present known.) Eigenstates of 6*4 are defined by 

S4:\a) = fj.a\a) Ma=^ae~*^°, (7) 

((a I a) ~ 1) where ^q, = ±1 is the eigenvalue of S4. Inserting complete sets of states in the second line 

of eq. (^) gives the third line where the amplitude Aap = {0\B\a){P\Bw\0), and /iat — exp(— A/at), 
Ha = — exp (—Ma). (We shall call A, a little artificially, the parity partner to the nucleon.) The time 
box length (T, here = 24) has been assumed to be so large that there are no finite time-size affects 
(although these can be easily taken into account). The parity partner arises because on the lattice 
inversion Is : Xi ^ —Xi is not parity (P) as it effects the fiavour degrees of freedom. (One can show that 
P = /sS4.) As a consequence hadron operators in general do not have a definite parity. 

As an example of a correlation function we show in Fig. nl a baryon correlation function. The typical 




Figure 1: A typical baryon correlation function plotted against t for (3 = 5.35, m = 0.01. Boxes denote 
negative values. The fit interval is from i = 2 to 22. 



zigzagging of the results (and fit) is caused by the parity partner. To pick out the ground states we must 
take t ^ 00 and/or choose an operator with a good overlap with this state. Experience has shown that a 
non-local operator (the 'wall' and hence the subscript W) has rather good overlap properties with the 
ground state for very moderate t, so this has been used as the source; the sink is taken as a local baryon 
operator. The mass obtained this way agrees with that which comes from using a local baryon operator 
for both sink and source. 

Measuring aTrN 

Practically there are several possibilities open to us for the evaluation of the matrix element. The 
easiest is simply to differentiate eq. (0). As0 = — XX'S'4 we have 

m^\p^{N\mxx\N) = a^N, (8) 

an application of the Feynman-Hellmann theorem. Thus all we need to do is to measure Mn for different 
masses m (at constant /3) and numerically estimate the gradient. In Fig. ^ we show such a procedure. 
(At present we have simply fitted a straight line through all points.) We find 

a^N - 11.9(8)mU=o.oi ~ 0.12(1), (9) 



the staggered formalism the operator expression for nixx has a more comphcated appearance than that given here, 



Figure 2: The Nucleon mass, Mjv, plotted against m for /3 = 5.35. 



which gives 



11.9TO 



M 



N 



11.9m + 0.64' 



0.15, 



(10) 



which is to be compared with the experimental result of u-^n /J^In ~ 0.06 ~ 0.05. The numerical result is 
much larger, but presumably this simply indicates that we have used much too large a quark mass in our 
simulation (which at present we, and everybody else, cannot avoid doing due to computer limitations). 

We now wish to estimate the valence and sea contributions. This is technically more complicated and 
involves the evaluation of a 3-point function, |L1| 



C{t;r) = {B{t)mxx{T)Bwm. 
This may be diagrammatically sketched as the sum of two terms: 



(11) 




(a) 



Bw \ j B (b) 

t 

Diagram (a) is the connected (or valence) piece; the nucleon quark lines are connected with the mass 
insertion operator. Diagram (b) represents the disconnected (or sea) piece. 

For the connected piece we have fixed t to be 8 or 9 and then evaluated C{t\ t) as a function of r. 
The appropriate fit function is 

C{t-T)^ A„/3(a|mxxl/3)ML">^ » t » r » 0. (12) 

The Aaf3 and are known from 2-point correlation functions. We see that when a — (3 — N we have 
the matrix element that we require: {N\mxx\N) ■ However there are other terms which complicate the 
fit: (A|mxx|A) and the cross terms (A|mxx|A^) = ^i^^'^ {N\mxx\^) ■ As can be seen from eq. ( p!2| ) 
these cross terms are responsible for oscillations in the result. To disentangle the wanted result from 
(A|mxx|A) we need to make a joint fit to the i = 8 and < = 9 results. At present we have not done this. 




Figure 3: The ratio of the 3-point connected correlation function to the 2-point correlation function 
against r for /? = 5.35, m = 0.01. The points used in the fit are joined together. (The dotted lines just 
extend the fit to neighbouring points.) 

but just checked that separate fits give consistent results. In Fig. |^ we show preliminary results for t ^ 8. 
Up until now we have only evaluated about a quarter of our available configurations, so the final results 
should have reduced error bars. We find 



« 0.08(2), 



(13) 



with about the same result for the A matrix element. The cross term is smaller, roughly 0.02. 

Finally we have attempted to estimate the disconnected term using a stochastic estimator, ||l^. This 
is costly in CPU time. One can improve the statistics by summing the 3-point correlation function over 
r; this is then equivalent to a differentiation of the 2-point function with respect to rn. This gives 



ir>t> 0. 



(14) 



a=N.A 



In Fig. ^ we show the disconnected part of the 3-point correlation function. As expected the quality of 




Figure 4: The ratio of the 3-point disconnected correlation function to the 2-point correlation function 
against t. 100 sets of Gaussian random numbers were used for the stochastic estimator. To improve the 
signal we have averaged over the points t and T — t. 



the data is poor - optimistically we see a slope. A simple linear fit gives 

0.05(4). 



(15) 



Discussion 

We see that (roughly at least) eq. (^) is consistent with eqs. (|l3|,^5|) and that 

^ « 1.5, (16) 

which is to be compared with the experimental value of about 2.2 ^ 1.8. Although we can draw no 
firm conclusions at present our result tentatively indicates that the valence part of the ttN sigma term is 
slightly larger than the sea term. 

Comparing our results with other work obtained using dynamical fermions, Gupta et. al, |Q who use 
2 flavours of Wilson fermions, find CTttA' /cjjjv ~ 2 ~ 3, which indicates a somewhat larger sea component 
in the ttN sigma term. On the other hand Patel, [Esl, using results from jlj] for 2 staggered flavours 
finds for the ratio 1.5 ~ 2.0, while Bernard et. al., [|5|, have ~ 2.0. These, like our result, seems to be 
lower than for the Wilson fermion case. 

We would also like to emphasise that although lattices can give a first principle calculation, at present 
one is not able to do this. Technically the fit formulae that we employ, eqs. (|l^, |lj) are true only for the 
complete correlation function. However we have used them for either the connected or the disconnected 
part separately. The best way to circumvent this problem is to make simulations at different strange 
quark masses; differentiation as in eq. (^ would then give directly ms{N\ss\N), the strange content of 
the nucleon. However this calculation is not feasible at the present time. This should be the ultimate 
goal of lattice simulations, as other recent (non-lattice) theoretical results, ||^, have hinted that perhaps 
the strange quark content of the nucleon is not as large as supposed, previous results being explained 
by a combination of factors, such as S 7^ aT^N and higher order corrections to first order perturbation 
theory. (Indeed there are already tantalising lattice indications that this may be so, |]l6[ 0.) 

In conclusion we would just like to say that lattice results at present are generally in qualitative agree- 
ment with other theoretical and experimental results. However much improvement in the calculations is 
required to be able to make quantitative predictions. 
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